UH-511-970-00 
hep-ph/0008160 



HIGGS TRANSVERSE MOMENTUM AT THE LHC 



C. BALAZS 

Department of Physics and Astronomy, University of Hawaii, 
Honolulu, HI 96822, U.S.A. 

Resummed QCD corrections to the transverse momentum (Qt) distribution of the Standard 
Model Higgs boson, produced at the CERN Large Hadron CoUider, are presented. The small 
Qt factorization formalism is reviewed, which is used to extend the standard hadronic fac- 
torization theorem to the low Qt region, with emphasis on the matching to the standard 
hadronic factorization. Comparison of the Qt predictions from the extended factorization 
and the parton shower method is performed. 



1 Introduction 

One of the fundamental open problems of the Standard Model (SM) is revealing the dynamics of 
the electroweak symmetry breaking (EWSB). The physical remnant of the spontaneous EWSB, 
the Higgs boson, is the primary object of search at present and future colliders.lll At the 14 TeV 
center of mass energy Large Hadron Collider (LHC) at CERN, the SM Higgs boson will mainly 
be produced in proton-proton collisions through the partonic subprocess gg (via top quark loop) 
HX. To experimentally understand the Higgs signature, to enhance the statistical significance 
of the signal, and to measure its basic properties (mass, lifetime, spin, charge, couplings), it is 
necessary to determine the transverse momentum {Qt) of the Higgs bosons. 

To reliably predict the Qt distribution of Higgs bosons at the LHC, especially for the low to 
medium Qt region where the bulk of the rate is, the effects of the multiple soft~gluon emission 
have to be included. This, and the need for the systematic inclusion of the higher order QCD 
corrections require the extension of the standard hadronic factorization theorem to the low Qt 
region. With a smooth matching to the usual factorization formalism, it is possible to obtain a 
prediction in the full Qt range. Since Monte Carlo event generators are heavily utilized in the 
extraction of the Higgs signal, it is crucial to establish the reliability of their predictions. Thus, 
results of the parton shower formalism are also compared to those of the analytic calculation. 



2 Low Qt Factorization 

In this section the low transverse momentum factorization formalism and its matching to the 
usual factorization is described. When calculating fixed order QCD corrections to the cross 
section dcj / dQdQ^dy for the inclusive process pp — > HX., the factorization theorem is invoked. 
The relevant variables Q, Qt, and y, are the invariant mass, transverse momentum, and rapidity 



of the Higgs boson, respectively. The standard factorization 



d(^=Y. J '^-^fji/hAxi,Q)daj,j^ fh/h2(x2,Q), (1) 

a convolution in the partonic momentum fractions xi and X2, fails when Qj- <^ Q as a result 
of multiple soft and soft+collinear emission of gluons from the initial state. The ratio of the 
two very different scales in the partonic cross section (Jj^j^ (identified by the partonic indices 
ji), produces large logarithms of the form ln((5^/Q|n) (being singular at Qt = 0), which are 
not absorbed by the parton distribution functions fj/h: unlike the ones originating from purely 
collinear parton emission. As a result, the Higgs Qt distribution calculated using the conven- 
tional hadronic factorization theorem is unphysical in the low Qt region. 

To resolve the problem, the differential cross section is split into a part which contains all 
the contribution from the logarithmic terms (W), and into a regular term (Y): 

dQdQ\dy " ^ ^^^^2). (2) 

Since Y does not contain potentially large logs, it can be calculated using the usual factorization. 
The W term has to be evaluated differently, keeping in mind that failure of the standard fac- 
torization occurs because ii neglects the transverse motion of the incoming partons in the hard 
scattering. As it is proven Hfl, small Qt factorization gives the cross section as a convolution of 
transverse momentum distributions 

W{Q,Qt,xx,X2) = J2 J d'^^T Cj^/h^{Q,kT,xi)Ej^j^{Q,QT)Cj^/h2{Q,QT -kT,X2). (3) 

Here Hj^jj is a hard scattering function, and Cj/h a^^e partonic density distributions depending 
on both longitudinal (x) and transverse (/ct) momenta, and on the scale of the factorization 
which is set equal to the hard scale Q. The convolution simplifies to a product in the Fourier 
conjugate, i.e. transverse position (6) space 

W{Q,b, xi,X2) =Cj^/h^{Q,b,xi)nj^j2{Q,b)Cj^/h2{Q,b,X2), (4) 

where W, Cj/^ and "Hjija the Fourier transforms of W, Cj^h a^d Hj^jj. The generalized parton 
distributions Cj/^, together with the hard scattering function, satisfy an evolution equation 
somewhat similar to the usual DGLAP equations. 

The evolution equation, for the production of a colorless boson, takes the formll 

^^W{Q,b,xi,X2) = - ^[A{as{fi),Ci) + B{as{Q),Ci,C2)], (5) 

where as is the strong coupling constant. It is customary to choose the renormalization scales 
arising in the evolution equation such that Ci = 2e~'^^ = Cq and C2 = 1. The solution of the 
above evolution equation leads to the expression! 

W{Q,b,x^,X2)=Cj^/hAQ^b,xi)e-^^'^''*^C^^/^,^{Q,b,X2), (6) 

with the Sudakov exponent 



(7) 



which resums the large logarithmic terms. The partonic recoil against soft gluons as well as 
the intrinsic partonic transverse momentum are included in the modified parton distributions 



^a/hix,b,Q). 



(8) 



The A and B functions, and the Wilson coefficients Cja are free of logs and safely calculable 
perturbatively as expansions in the strong coupling 



n=l 



(9) 



The process independent non-perturbative functions J^a/h^ describing long distance transverse 
physics, are extracted from low-energy experiments.! 

In the large Qt region, where Qt ^ Q, the standard factorization theorem is applicable. The 
matching of the small Qt region to the large Qt result is achieved via the Y piece. To correct 
the behavior of the resummed piece in the intermediate and high Qt regions, it is defined as the 
difference of the differential cross section calculated from the standard factorization formula at 
a fixed order n of perturbation theory and its Qt <C Q asymptote:0 



Y{Q,Qt, Xi,X2) 



dQ2 dQ^dy dQ"^ dQ^dy 



(10) 



Using this definition, the cross section to order is written as 



da 



dQ^ dQldy 



W{Q,Qt, Xi,X2) + 



dQ^ dQldy dQ^ dQldy 



(11) 



At low Qt-, when the logarithms are large, the asymptotic part dominates the Qt distribution, 
and the last two terms cancel in Eq. (0), leaving W well approximating the cross section. At 
Qt values comparable to Q the logarithms are small, and the expansion of the resummed term 
cancels the logarithmic terms up to higher orders in 05. Q In this situation the first and third 
terms nearly cancel and the cross section reduces to the fixed order perturbative result. After 
matching the resummed and fixed order cross sections in such a manner, it is expected that 
the normalization of the cross section ( |TT1) reproduces the fixed order total rate, since when 
expanded and integrated over Qt it deviates from the fixed order cross section a^"'^ only in 
higher order terms. Further details of the low Qt factorization formalism and its application to 
Higgs production can be found in the recent literatureB 



3 Higgs Qt at the LHC 

The low Qt factorization formalism, described in the previous section, is utilized to calculate 
the QCD corrections to the production of Higgs bosons at the LHC. In the low Qt region 
this calculation takes into account the effects of the multiple-soft gluon emission including the 
Sudakov exponent S and the non-perturbative contributions J-a/h- In Sudakov exponent 
the A^'^\ A^'^\ and B^^^i coefficients are included. The 0{a^g) virtual corrections are also taken 
into account by including the Wilson coefficient C^gg , which ensures the 0{a\) total rate. By 

"To prevent evalu ation of the Su dakov exponent in the non-perturbative region, the impact parameter 6 was 
replaced by b, = b/ ^1 + (fe/femax)^. m 
''The expression of the Y term for Higgs production can be found elsewherejj 
'^The cancellation is higher order than the order at which the singular pieces were calculated. 
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Figure 1: Higgs boson transverse momentum distributions calculated by ResBos (curves) and PYTHIA (his- 
tograms). The default (middle) ResBos curve was calculated with the canonical choice of the renormalization 
constants, and the other two with doubled (lower curve) and halved (upper curve) values of Ci and C2. For 
PYTHIA, the original output with default input parameters (dashed), the same rescaled by a factor of K = 2 
(dash-dotted), and a curve calculated by the altered input parameter value Qmax = s (dotted) are shown. The 
lower portion, with a logarithmic scale, also shows the high Qt region. 

matching to the 0{ag) fixed order distributions a prediction is obtained for the Higgs production 
cross section in the fuU Qt range which is vahd up to 0(q|.). The details of this calculation 
are given in an earlier work. @ The analytic results are coded in the ResBos Monte Carlo event 
generator. Q 

Fig. 1^ compares the Higgs boson transverse momentum distributions calculated by ResBos 
(curves) and by PYTHIAeI (histo grams from version 6.122). The middle solid curve is calculated 
using the canonical choice for the renormalization constants in the Sudakov exponent: Ci = Co, 
and C2 = 1. To estimate the size of the uncalculated B^"^^ term, these renormalization constants 
are varied. The upper solid curve shows the result for Ci = Co/2,C2 = 1/2, and the lower 
solid curve for Ci = 2Co,C2 = 2. The band between these two curves gives the order of the 
uncertainty originating from the exclusion of B^'^\ The typical size of this uncertainty, e.g. 
around the peak region, is in the order of ±10 percent. The corresponding uncertainty in the 
total cross section is also in the same order. This uncertainty is larger than the uncertainty 
arising from the non-perturbative sector of the formalism (which was estimated to be less than 
5 percent in the relevant Qt region). 

4 Comparison to Parton Showers 

Multiple soft-gluon radiation from the initial state can also be treated by the parton shower 
technique.u This approach is based on the usual factorization theorem, giving the probability 
V = e~^^^^ of the evolution from the scale Qo to Q, with no resolvable branchings by the 
exponent 

5(03)=^ *^5M/'<i.P„_a.), (12) 

Jq2 2tt Jo 
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which is defined in terms of the DGLAP sphtting kernels Pa^bc{z)- The formalism can be 
extended to soft emissions as well by using angular ordering. The distinct difference between the 
Sudakov exponents of the low Qt factorization and the parton showering approach is apparent. 
A comparison can be made based on the qualitative argument that parton showering resums 
leading logs which depend only on the given initial state. These logs correspond to the logs 
weighted by the A function in Eq.(|^. 

Fig. H shows that the shape of the PYTHIA histogram agrees reasonably with the resummed 
curve in the low and intermediate Qt (~ 125 GeV) region. For large Qt, the PYTHIA prediction 
falls under the ResBos curve, since ResBos mostly uses the exact fixed order 0(a|) matrix ele- 
ments in that region, while PYTHIA still relies on the multi-parton radiation ansatz. PYTHIA 
can be tuned to agree with ResBos in the high Qt region, by changing the maximal virtuality 
a parton can acquire in the course of the shower (dotted curve). In that case, however, the low 
Qt region will have disagreement. 

Since showering is attached to a process after the hard scattering takes place, and the parton 
shower occurs with unit probability, it does not change the total cross section for Higgs boson 
production given by the hard scattering rate. Thus, the total rate is given by PYTHIA at 
0{a'g). In Fig. |l] the dashed PYTHIA histogram is plotted without altering its output. For 
easier comparison the default PYTHIA histogram is also plotted after the rate is multiplied by 
the factor K = 2 (dash-dotted). A detailed comparison of the results for Higgs boson production 
from ResBos and from event generators based on the parton shower algorithm is the subject of 
a separate work.EI 
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